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Abstract 

We investigate the possibility of observing the direct CP violation in 
the decay modes B~ — > D°D~ and D°D~ within the Standard Model. 
Including the contributions arising from the tree, annihilation, QCD 
as well as electroweak penguins with both time- and space-like compo- 
nents, we find that the direct CP asymmetry in B~~ — > D°D~ is very 
small ~ 0.2 % but in B~ — > D°D~ decay it can be as large as 4%. 
Approximately 10 7 charged B mesons are required to experimentally 
observe the CP asymmetry parameter for the later case. Since this is 
easily accessible with the currently running B factories, the decay mode 
B~ —* D°D~ may be pursued to look for CP violation. 

PACS Numbers. : 11.30 Er, 13.25 Hw 



1 Introduction 

CP violation is one of the least understood phenomena in particle physics 
|], H], [|, although it was observed in K° — K° mixing system more than 35 
years ago. In the standard model (SM), CP violation arises from a complex 
phase in the Cabibbo-Kobayashi-Maskawa quark mixing matrix 0. Outside 
the Kaon system, decays of B mesons provide rich ground for investigating CP 
violation @, §. Within the SM, the CP violat ion is often characterized by the 
so-called unitarity triangle JR. By measuring CP violating rate asymmetries in 
B decays, one can extract a, (3 and 7, the three interior angles of the unitarity 
triangle. The sum of these three angles must be equal to 180° in the SM with 
three generations. At present we are at the beginning of the 5-factory era 
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in particle physics, which will provide us valuable insights to understand the 
phenomena of CP violation. One of the main programmes of the presently 
running and the upcoming B factories is to measure the size of CP violation 
in as many B decay modes as possible so as to establish the pattern of CP 
violation in various B decays. Among the most interesing B decay channels, 
the "gold plated" mode B d — > J/ipK a , allows the determination of the 
angle (3 of the unitarity triangle of CKM matrix. Recent measurement of CP 
asymmetry in the B° — > J/ipK° and other related processes e.g. ip'K , r] c K° 
etc by the BELLE gj and BaBar detectors at the KEK and SLAC B 



factories together with the earlier measurement of CDF |TlJ| constitute the 
first significant signal of CP violation outside the neutral Kaon system. 

While the most promising proposal for observing CP violation in the B- 
system involves the mixing between neutral B mesons |l| , the decays of charged 
B mesons are also of particular importance for establishing the detailed nature 
of CP violation. Since charged B mesons cannot mix, a measurement of the 
CP violating observable in these decays would be a clear sign of "direct CP vi- 
olation" which has been searched for in i^-system for quite long with indefinite 
success. Only recently, such kind of CP violating effect has been observed in 
the K system by the NA48 [|12] and KTeV |13j collaborations. For the bottom 
meson case usually the charmless rare B decay modes are preferred to study the 
direct CP violation as these decay modes proceed with more than one Feynman 
diagrams. In this paper we would like to look for some additional decay chan- 
nels which could help us in establishing the presence of CP violation as quickly 
as possible. For this purpose we investigate the direct CP violating effecs in 
the decays of charged B mesons to two charmed mesons i.e. B~ — > D°D~ and 
D°D~. It is worth emphasizing that these decay modes are flavor self tagging 
processes which should be favored for experimental reconstuctions. The decay 
mode B~ — > D°D~ has already been observed experimentally with a branch- 
ing ratio (1.3 ±0.4)% and the upper limit for B~ — > D°D~ channel is found to 



be < 6.7 x 10 3 [1_4|. These decay modes which are described by the the quark 



level transitions as b — > ccq (q = s/d for D~ / D~ in the final state) proceed 
through three distinct type of flavor topologies. These are : the color allowed 
but Cabibbo suppressed tree, annihilation and the QCD as well as electroweak 
penguin diagrams. To get significant direct CP violation one would require 
two interferring amplitudes of comparable strengths, with different strong and 
weak phases. The weak phases arise from the superposition of various penguin 
contributions and the usual tree diagrams. The strong phases are generated by 
the perturbative penguin loops (hard final state interaction) [OJ or final state 
interactions involving two different isospins. Since the decay modes we consid- 
ered here are single isospin channels i.e. the final states D°D~ and D°D~ are 
with isospin 1 = 1/2 and 1 respectively the second type of FSI strong phase 
differences are absent for these channels. Therefore at the first sight it appears 
that direct CP violating effects in these channels would be negligibly small as 
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the tree contribution dominates over the other diagrams and thus have been 
overlooked in the literature. But detailed calculation shows that it is indeed 
not so. In fact the CP violating effects in B — > D°D~ channel can be as large 
as few percent level which can be experimentally accessible in the first round 
of B factories. The reason for the existence of such a significant CP violating 
parameter may be due to the fact that the tree diagram for b — > ccd transi- 
tion is although colour allowed, it is doubly Cabibbo suppressed, and hence 
its magnitude is not very much larger than the penguin contributions. CP 
violating effects in the decays of neutral B meson into double charmed mesons 
have been extensively studied in Refs. [16-19], where it has been shown that, 
these channels can be used as an alternative method to the J/ipK s mode for 
the extraction of the angle j3. 

In our analysis, we use the standard theoretical framework to study the 
nonleptonic B~ — > D°D~(D~) decay modes, which is based on the effective 
Hamiltonian approach in conjuction with the factorization hypothesis. The 
short distance QCD corrected Hamiltonian is calculated to next-to-leading 
order. The renormalization scheme and scale problems with factorization ap- 
proach for matrix elements can be circumvented by employing the scale and 
scheme independent effective Wilson coefficients. In the literature the con- 
tributions of space-like penguins are neglected assuming form factor suppres- 
sion. But as pointed out in Ref. |20| the effect of space like penguin ampli- 
tudes can be remarkably enhanced by the hadronic matrix elements involving 
(V — A)(V + A) or (S + P)(S — P) currents. Therefore we have included the 
space and time like contributions of both QCD and EW penguins, the anni- 
hilation contribution in addition to the dominant tree diagrams. Assuming 
the factorization approximation, the matrix elements of the tree and time-like 
penguin diagrams have been calculated in the BSW model JH], whereas for 
the evaluation of the matrix elements of the space and annihilation diagrams 
we have employed the Lepage and Brodsky model p2[ . 

The paper is organized as follows. In section II we briefly discuss the ef- 
fective Hamiltonian together with the quark level matrix elements and the nu- 
merical value of the Wilson coefficients in the effective Hamiltonian approach. 
Assuming the factorization approximation, the matrix elements of tree and 
time-like penguins are evaluated in the BSW model and for the space-like and 
annihilation diagrams we use LB (Lepage and Brodsky) model. Determina- 
tion of the CP violating asymmetry is presented in section III and section IV 
contains our conclusion. 



2 Framework 

The effective Hamiltonian 'Heff for the decay modes B~ — > D°D~ and D°D~ 
which are described by the quark level transitions b —>■ ccq (where q = s for 
the former and d for the later) have three classes of flavour topologies : the 
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dominant tree, annihilation, and both QCD as well as electroweak penguins 
given by |§ 

n eff = ^{k[cMo u M + cMo u m] + \ e [c 1 fa)o°M + c 2 (//)o 2 »] 

10 . 

+ (A u + A c )^ Ci (/i)Oi(/i) +h.c. , (1) 

i=3 J 

where X u = V u bV* q and A c = V c bV* q and Ci(fi) are the Wilson coefficients 
evaluated at the renormalization scale /i. The four fermion operators Oi_io 
are given as 



O" = (ub) V -A(qu)v-A , 0% = (u a bf3) V - A (qi3U a )v-A , 

0\ = (cb) v ^ A (qc)v- A , O2 = (c a b[3) V - A (qpc a ) V - A , 

03(5) = {qb) v _ A Y^{q' q')v-A{v+A) , 



04(6) = (qMv-Aj2(^'a] 



V-A(V+A) , 
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7 (9) = ^{qb)v-A^Ze q >{q q') v+A (v-A) , 
g' 

3 

8 (io) = 2(^« & / 3 )^-A^Z e 9'(^^)v+A(v-A) , (2) 

q' 

where Oi i2 are the tree level current- current operators, the QCD pen- 

guin operators and O-j-iq the EW penguin operators. (qiq2){v±A) denote the 
usual (V ± A) currents. The sum over q' runs over the quark fields that 
are active at the scale /i = 0{m b ) i.e. (q f G u,d, s,c,b). The Wilson coef- 
ficients depend (in general) on the renormalization scheme and the scale [i 
at which they are evaluated. In the next to leading order their values ob- 
tained in the naive dimensional regularization (NDR) scheme at fi = rrib(mb) 
as H d = 1.082, c 2 = -0.185, c 3 = 0.014, c 4 = -0.035, c 5 = 0.009, c 6 = 
-0.041, c 7/a = -0.002, c 8/a = 0.054, c 9/a = -1.292 and c w/a = 0.263. 

However, the physical matrix elements (PiP 2 |7i e / f\B) are obviously inde- 
pendent of both the scheme and the scale. Hence the dependence on the 
Wilson coefficients must be compensated by a comensurate calculation of the 
hadronic matrix elements in a nonperturbative framework such as lattice QCD. 
Presently, this is not a viable strategy as the calculation of the matrix elements 
{PiP2\Oi\B) is beyond the scope of the current lattice technology. However, 
perturbation theory comes to (partial) rescue; with the help of which one-loop 
matrix elements can be rewritten in terms of the operators and the effective 
Wilson coefficients c/ which are scheme and scale independent : 

(qq'q'\H eff \b) = £ cf f {^) (qq> q' \0 3 \b)^ . (3) 
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The effective Wilson coefficients c% (p) may be expressed as [2| 



r eff 
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(4) 



where f T and 7 < -°^ T , the transpose of the matrices f and arise from the 
vertex corrections to the operators 0\ — 0\q derived in which are explicitly 
given in Ref. |26| 

The quantities C t , C p , C e and C g are arising from the peguin type diagrams 
of the operators Oi,2, the QCD penguin type diagrams of the operators 3 —0 6 , 
the electroweak penguin type diagrams of 0\ i and the tree level diagrams of 
the dipole operator O g respectively, which are given in the NDR scheme (after 
MS renormalization) by 



c t = - (t^KJ + Y t ^ mc) ) Cl 

C p = [G(m s ) + G(m b )}c 3 + J2 G(m<) (e* + c 6 ) 



C„ 



2rrib 



i=u,d,s,c,b 

?/ f , c e / f = -1.043 
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C e = ~ (yG(m u ) + yG(m c )\ ( Cl + 3c 2 ) 



2 

G(m q ) = -- G{m q ,k,p) (5) 



G(m, fc, //) = — 4 / <ir x(l — x) In , (6) 

Jo \ p 2 J 

It should be noted that the quantities C t , C p C e and C g depend on the 
CKM matrix elements, the quark masses, the scale p and k 2 , the momentum 
transferred by the virtual particles apearing in the penguin diagrams. In the 
factorization approximation there is no model independent way to keep track 
of the k 2 dependence; the actual value of k 2 is model dependent. From simple 
kinematics |27|] one expects k 2 to be typically in the range 



m b 2 . m 



, <r<-±. (7) 

Since the branching ratio and the CP asymmetry depend crucially on the 
parameter k 2 , here we would like to take a specific value for it based on valence 
quark approximation instead of the conventionally used value k 2 = m|/2. As 
discussed in Ref. |HJ the averaged value of the squared momentum transfer 
for B~(bu) — > D°(cu)D~ (qc) is given as 

(k 2 ) =m 2 + m 2 q - 2m b E q (8) 
where the energy of the quark q in the final D~ particle is determinable from 



E q + ^JE 2 -m 2 q + m 2 + ^J4E 2 - Am 2 q + m 2 c = m b (9) 
for time-like penguin channels; or from 



E q + ^JE 2 -m 2 + ml = m b + m u (10) 

for space-like penguin diagrams. m&, m q and m c denote the masses of the 
decaying 6-quark, daughter g-quark and the c-quark (created as cc pair from 
the virtual gluon, photon or Z-particle in the penguin loop). For numerical 
calculations, we have taken the CKM matrix elements expressed in terms of the 
Wolfenstein parameters with values A = 0.815, A = sin 6 C =0. 2205, p = 0.175 
and 7] = 0.37 PBf. The choice of p and rj correspond to the CKM triangle : a = 



91°, (3 = 24° and 7 = 65°. At scale p ~ m b , we use the current quark masses 
as m u {m b ) = 3.2 MeV, m d {m b ) = 6.4 MeV, m s {m b ) = 90 MeV, m c {m b ) = 
0.95 GeV and m b (m b )= 4.34 GeV. With the specific value of k 2 obtained from 
Eqns. (8-10), we obtain the values of the effective renormalization scheme and 
scale independent Wilson coefficients for b —> s and b d transitions as given 
in Table- 1. 
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Now we want to calculate the matrix element (D~ D°\Oi\B~) using the fac- 
torization approximation, where Oi are the four quark current operators listed 
above. In this approximation, the hadronic matrix elements of the four quark 
operators (cb)(y-A)(Qc)(v-A) split into the product of two matrix elements, 
(D°\(cb)(y-A)\B~) and (D~\(qc)(y-A)\fy where Fierz transformation has been 
used so that flavor quantum numbers of the currents match with those of 
the hadrons. Since Fierz rearranging yields operators which are in the color 
singlet-singlet and octet-octet forms, this procedure results, in general, in ma- 
trix elements which have the right flavor quantum numbers but involve both 
singlet-singlet and octet-octet current operators. However, there is no exper- 
imental information available for the octet-octet part. So in the factorization 
approximation, one discards the color octet-octet piece and compensates this 
by treating N c , the numbers of colors as a free parameter, and its value is 
extracted from the data of two body nonleptonic decays. 

The matrix elements of the (V — A)(V + A) operators i.e. (Oq & 8 ) can be 
transformed into (V — A){V — A) form by using Fierz ordering and the Dirac 
equation, which are given as 

(D;D \O 6 \B-) = R q {D- q D°\0,\B-) (11) 

with 

2m D - 

R q = 7 \7~^ 1 7 > ( 12 ) 

(m b - m c )(m q + m c ) 
where the quark masses are the current quark masses. The same relation works 
for 8 . 

Hence, one obtains the transition amplitude for B~ — > D~D° and D~D° 
as (where the factor Gf/V% is suppressed) 

A(B~ - DJD°) = A M { (a A + a 10 + (a 6 + a 8 )R^X^ BD °' D ^ 

+ (a! + a 4 + an, + K + a 8 )R'^jX {B ' D ° D ^ 

+ A c | (a! + a 4 + ai + K + a R )R 8 ^ X( BD °> 

+ ( K a 4 + a 10 + (a 6 + a 8 )R^X^ D ° D ^ (13) 

A(B~ - J D' J D°) = A u |(a 4 + a 10 + (a 6 + a 8 ) J R (i )x( i?D0 ' D ^ 

+ [a, + a 4 + aio + (a 6 + a 8 )R' d ^X^ D ° D ~^ 

+ A c | (en + a 4 + aio + K + a 8 )R d ^j X (BD °' D_) 

+ (a 4 + a 10 + (a 6 + a 8 ) J R^X^ D °^)} (14) 



,D7) 
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where 



X ( BD °^ = (D;\(qc)\0)(D°\(cb)\B) 

x (b,d-do } = {D ° D - miQ) ( |(fi6)|S> . (15) 

X(BD°,D q ) d eno t es matrix elements of the tree and time like penguins where 
as X( B,D ° D i ) stand for the annihilation and space-like amplitudes. 

2m% 

R q ' = 7 T7 ; r , (16) 

(m q - m u )(m b + m u ) 

arises from the transformation of (V — A)(V+A) operators into (V — A)(V — A) 
form for space-like penguins. It should be noted that \ u = V u bV* s for B~ — > 
D°D~ whereas X u = V u bV* d for B~ — > D°D~ and similar expressions for A c . 

The coefficients at, ai ■ ■ ■ aio are combinations of the effective Wilson 
coefficients given as 

a 2i -i = + -^ff c t S «2i = c 2 { f + -^rjciU i = 1, 2 • • • 5 , (17) 

JMc JMc 



where N^" is the effective number of colors treated as free parameter in order 
to model the nonfactorizable contributions to the matrix elements and its value 
can be extracted from the two body nonleptonic B decays. A recent analysis 
of B — > Dtt data gives N^" ~ 2 |28||. Therefore, in our analysis, we take two 



c 

sets of values for N^" i.e., N%** = 2 and N^'* = 3, which characterizes naive 
factorization. 

The factorized hadronic matrix elements are evaluated using the BSW 
model [21 1, which are given as 

XiBD^u-) = ifDqF B D (ml g )(ml - ml ) (18) 

The matrix element of the annihilation and space-like penguins are given 



as 



(D°D-\(qu)(ub)\B-) =zf B fl(m 2 B ) 



2 2 m D q ~ m D 2 

m D q ~ m D° ~ - — - m B 



(19) 



where the value of the annihilation form factor is given as f1{m 2 B ) = i\ l o'na s j B lrr? B 
I 

After obtaining the transition amplitude, the branching ratio is given as 

Ipl \A(B- ^ D°D-)\ 2 

BR=^\ 1 -^ 2iL, (20) 

8irm 2 B r v ; 

where |p| is the momentum of the emitted particles and T is the total decay 
width. 
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Using eqns (|TB|)-(0) we obtain the transition amplitude (in the unit of 

G F /y/2) as 



A(B- 



D°D;) = A u (0.1898 - z0.6483) + A c (0.1889 + i4.418) 
[A u (0.20i9 - z0.6817) + A c (0.201 + i4.698)] , (21) 



A{B~ -> D°D~) = A u (0.2259 - z0.5616) + A c (0.2259 + z4.8185) 
[A M (0.2393 - i0.589) + A c (0.2393 + z'5.124)] , 



:22) 



where we have used the decay constants (in MeV) as fu s = 280, fu = 300 
I4| and /s=180 In the above equations, the upper values correspond to 
N^ff = 2 and the lower bracketed values to N^* = 3. 



3 CP Violating Asymmetry 

For charge B T decays, the CP violating rate asymmetries in partial decay rates 
are defined as 

_ T(B- -> /-) - T(B+ -> /+) 
cp ~ T(B- -> /-) + T(B+ -> /+) • 1 J 

As these decays are all self tagging the measurement of these CP violating 
asymmetry is essentially a counting experiment in well defined final states. 
Their rate asymmetries require both weak and strong phase differences in in- 
terfereing amplitudes. The weak phase difference arises from the superposition 
of amplitudes from various tree (current-current) and penguin diagrams. The 
strong phase which are needed to obtain nonzero values for generated 
by absorptive parts in penguin diagrams (hard final state interactions). 

For the B meson decaying to a final state / and the charge conjugated 
B~ — > / we may, without any loss of generality, write the transition amplitude 
as 

A(f) = \ u A u e i5 « + \ c A c e i5c (24) 

A(f) = \* a A u e iSu + X*A c e i5c (25) 

where Aj = VibV*, A u and A c denote the contribution from penguin operators 
proportional to the product of CKM matrix elements A u and A c respectively. 
The corresponding strong phases are denoted by S u and S c respectively. 
Thus the direct CP violating asymmetry is given as, 

-2 Im(A u A*) lm{A u A* c ) 

ct, 



" " \KAu\ 2 + \X C A C \ 2 + 2 Re(A u A*) Re(A u A 
2 sin 7 sin(5 M — 8 C ) 



^| + |^|+2cos 7 cos(5 u -5 c ) ' (26) 
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where the weak phases entering in the b — > s/d transition is equal to (—7), as 
we are using Wolfenstein approximation in which A c has no weak phase and 
the phase of \ u is —7. The strong phase (5 U — S c ) is caused by the final state 
interactions. 

The strong phases are given by, 

sin(5 u - S c ) = (lmA u ReA c - lmA c ReA u ) (27) 

cos(5 u - S c ) = 1 {ReA u ReA c + lmA u lmA c ) (28) 



4 Conclusion 

Using the next-to-leading order QCD corrected effective Hamiltonian, the scale 
and scheme independent Wilson coefficients, we have systematically studied 
the two charm hadronic decay modes B~ — > D°D~ and D°D~ within the 
framework of generalized factorization. The nonfactorizable contributions are 
parametrized in terms of N*" , the effective number of colors. For numerical 
calculations, we have used two different sets of values for this parameters: (i) 
N e Jf = 2, (ii) N e Jf = 3 which holds for naive factorization. The existence 
of a direct CP violating rate asymmetry requires two interferring amplitudes 
having different CP nonconserving weak phases and CP conserving strong 
phases. The former may arise either from the Standard Model CKM matrix 
or from new physics while the latter may arise from the absorptive part of a 
penguin diagram or from final state interaction effects of two different isospins. 
Since the channels we considered here are single isospin channels, the second 
class of strong phase differences do not arise for these channels. In our analysis, 
the weak phases are due to the CKM matrix and the strong phase differences 
arise due to absorptive part of penguin diagrams. The branching ratio and the 
CP violating asymmetry parameter are estimated using Eqs.(^) and (|26| ) and 
are presented in Table-2 

From the results we have observed the following: 

1. The predicted branching ratio for the decay mode B~ — > D°D~ agrees 
very well with the experimental value for Nf:** = 2, and the CP violating 
parameter for this mode is quite small. 

2. The branching ratio for the decay mode B~ — > D°D~ lies below the 
present experimental upper bound and the CP violating parameter for this 
mode is quite significant. The number of charged B mesons required to observe 
this CP violating signal to three standard deviation is given as Nb = 9 /(BR x 
a cp) ~ 7.9 x 10 6 , which is easily accessible with the currently running B 
factories. 

It has been emphasized in Refs. [16-19] that the neutral B meson decay 
modes to two charmed mesons can be used to measure the unitarity angle (3 
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as alternative to the gold plated mode B — > J/ipK. We argue further here 
that the mode B~ — ► D°D~ can be used to quickly settledown the search for 
observing direct CP violation outside Kaon system, if SM description of CP 
violation is correct or else could provide us a clear indication of the presence of 
new physics. It should be noted here that the decay mode is flavor self-tagging 
and hence experimentally favourable. Furthermore, since the branching ratio 
is O(10 -3 ) in our case, which is larger than the other competitive mode like 



B — > ttK where the branching ratio is O(10 ) [14] and the direct CP violation 
parameter a cp is significantly large i.e., we have obtained a cp around 4 % as to 
that of —1.4 % in ref JHJ, this decay mode is certainly a better candidate to 
observe CP violation in the first round of B-factory experiments. 

To summarize, since the modes we consider are direct decays and not time 
dependent, they may be observed in any experimental setting where large num- 
ber of B mesons are produced. Apart from the SLAC and KEK asymmetric B 
factories these include CLEO and hadronic B experiments such as HERA-b, 
BTeV, Collider Detectors at Fermilab (CDF), DO and CERN LHC-b or high 
luminosity Z factory. As we have used generalised factorization approxima- 
tion along with BSW model and Lepage and Brodsky model for penguin and 
annihilation contributions, we might have introduced certain uncertainties. 
Nevertheless, since the branching ratio obtained for B~ — > D°D~ matches 
very well with the experimental value, we point out here that the decay mode 
B~ — > D°D~ may also be pursued further in the first round of B-factory ex- 
periments (where it can easily be accessible) to observe direct CP violation or 
to provide us a hint for the presence of new physics. 
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Table 2: Branching ratio and CP Asymmetry in % for B — > D°D~,D°D~ 
decay modes. 



Decay modes 


Branching 
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CP 


Asymmetry 




iV c e ^' = 2 


iV c e// = 3 
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N%" = 2 


Nl" = 3 


B~ -> D°Dj 


1.29 x 1(T 2 


1.46 x 1(T 2 


(1.3 ±0.4) x lO- 2 
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0.18 


B~ -> D°D~ 
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9.86 x 1(T 4 
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